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Science Frontier Panels: 
^ ' Galaxies Across Cosmic Time ( GCT) 

Projects/Programs Emphasized: 
1. The Energetic X-ray Imaging Survey Telescope (EXIST); 
http : / / exi st . gsf c . nasa . gov| 
2. The Wide-Field X-ray Telescope (WFXT);'http : / /wf xt . pha . jhu . edul 
3. Panoramic Survey Telescope & Rapid Response System (Pan-STARRS); 
|http : //pan- S tar rs . if a . hawaii . edu/public/ 
4. The Large Synoptic Survey Telescope (LSST); [http : //Isst . org 



5. The Synoptic All-Sky Infrared Survey (SASIR); http : / /sasir . org 



Key Questions: 

1. What is the assembly history of massive black holes in the uni- 
verse? 

2. Is there a population of intermediate mass black holes that are 

the primordial seeds of supermassive black holes? 

3. How can we increase our understanding of the physics of accre- 
tion onto black holes? 

4. Can we localize sources of gravitational waves from the de- 
tection of tidal disruption events around massive black holes and 
recoiling binary black hole mergers? 



1 Introduction 

Dynamical studies of nearby galaxies suggest that most if not all galaxies with a bulge component 
host a central supermassive black hole (SMBH), and that the bulge and BH masses are tightly 
correlated dUSSI^IliQ]. This is referred to as the Mbh— cr* relation, where the velocity dispersion 
(a^) of bulge stars is a proxy of the bulge mass. Why the BH mass, which is concentrated on 
such small scales, should correlate with the bulge mass on such large scales is a deep puzzle. 
In the standard Lambda cold dark matter (ACDM) scenario, bulges supposedly grow by mergers 
whereas BHs grow from accretion of stars and gas (and perhaps dark matter [7]). Mergers induce 
gravitational perturbations that allow gas to rapidly infall fueling BH growth and ensuring at 
some level concurrent growth on large and small scales [9]. In what appears to be an efficient 
feedback mechanism [10], if the MBH lights up as an active galactic nucleus (AGN), then the 
growth of the BH can be halted as energy is pumped back into the gas reservoir quenches further 
growth [lH [12, [isi 14, 15, Even in the absence of r nerger- induced growth, random stellar 



encounters should serve to continually feed SMBH growth [llTl llSn. As a third channel for growth, 
BHs may grow via mergers with other BHs. As summarized by Treu et al. [19], in the absence 
of feedback, the slope of the Mbh — cr* relation could evolve in either sense depending on the 
respective growth histories. 

Clearly, gaining an unbiased (or at least a differently biased) understanding of the demographics 
[|2oll of BHs in galactic centers as a function of cosmic time is crucial for understanding the 
growth of structure in the universe. Yet the direct detection of SMBHs in the centers of quiescent 
galaxies via stellar and gas dynamics is severely limited by our ability to spatially resolve the 
sphere of influence of the central BH, Rsp^ ~ GMbh/ci'*- Thereby direct BH mass measurements 
are largely restricted to a limited number of nearby non-active systems for which accurate BH 
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masses can be obtained from dynamical measurements [12 ll. |22| . |23L l24l 12511 and active galaxies with 
intensive reverberation mapping campaigns 11261 1271. 12811 . 



2 Tidal Disruption Events: a tool for probing black holes 

Stellar dynamical models predict that once every 10^ — 10^ yr a star in a galaxy will pass within 
the tidal disruption radius of the central BH, < Rt = -R^(Mbh/^*)^^^, and will be torn apart 
by tidal forces [|29L [SOfl. After the star is disrupted, at least half of the debris is ejected from the 
system, and some fraction remains bound to the BH and is accreted. The fallback of debris onto 
the BH produces a luminous electromagnetic flare that is expected to peak in the UV/X-rays and 
radiate close to the Eddington luminosity 113 ll |32| . 13311 . 

The detection of a tidal disruption event (TDE) is unambiguous evidence for the presence of a 
central BH, and enables the detection of quiescent SMBHs - those not accreting, or at such low 
rates as to be undetectable. TDEs also provide a cosmic laboratory for studying the detailed physics 
of accretion onto BHs. The properties of the flare from the accreting stellar debris (luminosity, light 
curve, and spectral energy distribution) are dependent on the mass and spin of the BH. Thus, de- 
tailed observations of tidal disruption flares can provide an independent means of measuring 
the masses and spins of dormant BHs in distant galaxies. The disruption of compact objects, 
specifically white dwarfs, is also an interesting possiblity since it can only occurr around BHs with 
masses of order lO^M© or less and should have a distinct spectroscopic signature [34]. However, 
in the case of the most massive BHs, the tidal disruption radius is smaller than the event horizon 



and stars are swallowed whole without disruption (for a solar-type star Merit ~ 10 M©, n35\\ ). 



2.1 Progress to Date 

Although the volumetric rate of TDEs is low (10~^ — 10~^ yr~^ Mpc~^), about a dozen candidate 
TDEs have emerged in all-sky X-ray surveys (ROSAT and XMM-Newton Slew Survey) and the 
GALEX Deep Imaging Survey in the UV. 

The ROSAT AU-Sky Survey detected large amplitude (> 100), luminous (lO''^ - 10^^ ergs s^^) 
soft X-ray outbursts from 5 galaxies [36] that showed no previous evidence of AGN activity, and 
detection rates consistent with theoretical predictions for the rates of TDEs [37]. More recently 
the XMM-Newton Slew Survey BSL 13911 detected two galaxies with large amplitude (80 — 90), 



luminous (10 — 10 ergs s^ ) soft X-ray outbursts when compared to their ROSAT upper limits. 
Although two of these candidates (NGC 5905 and NGC 3599) were found to have optical spectra 
with narrow emission lines indicative of a low-luminosity AGN 11391. 14011 . the large amplitudes of 
their X-ray flares were more extreme than any previously known cases of AGN variability. These 
luminous, soft X-ray outbursts were best explained as the result of a sudden increase in fuel supply 
to an otherwise quiescent central BH due to the tidal disruption and accretion of a star. 
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The X-ray spectra of the candidate TDEs detected by ROSAT and XMM-Newton can be fitted 
with blackbody spectral energy distributions (SEDs) with temperatures of (6 — 12) x 10^ K, or 
extremely soft power-law spectra with F = 3 — 5, where /e oc E^^. The Rayleigh- Jeans tail of 
this blackbody emission should be prominent in the UV, as well as benefit from a strong negative 
K correction. In a systematic search of ~ 3 deg^ of the GALEX Deep Imaging Survey in the 
FUV and NUV two luminous (> 10^^ ergs s^^) flares were detected from galaxies at z = 0.37 and 
z = 0.33 with n o signs of AGN activity in their nuclei in the form of optical emission lines or hard 
X-ray emission (41, 421. 

The simultaneous optical difference imaging light curves of the UV flares, extracted from the 
CFHT Legacy Survey with an excellent cadence of days, are well described by the t^^/'^ power- 
law decay predicted both analytically and from numerical simulations to be t he p ower-law decline 
of the mass fallback rate of debris from a tidally disrupted star 1331 |43L l44l 1451] . The UV/optical 
SEDs of the flares are fitted with blackbody temperatures of ~ 5 x 10^ K, or power-laws with a = 
+1.0 — 1-1.5, where f,, oc z/". Soft X-ray emission was also detected during the flares, consistent 
with a second blackbody component with Tbb = (2 - 5) x lO^K. In Gezari et al. (2008) 142], they 
run a Monte Carlo simulation of their survey efficiency and find that the observed detection rate of 
the UV/optical candidates is consistent with the black-hole-mass dependent TDE rates from stellar 
dynamical models from Wang & Merritt (2004) [30] if the assumption is made that the fraction of 
flares that radiate at the Eddington luminosity decreases with increasing BH mass 11321] . and that 
stars are spun-up to near break up when disrupted 11461] . 



3 Future Prospects 

The detection of the light curves of tidal disruption events for the first time in the optical has 
encouraging implications for the next wave of high-cadence, optical/infrared time-domain surveys, 
which will be covering a large enough volume to detect tens of tidal disruption events per year. The 
detailed optical light curves and colors measured from high cadence difference imaging will be 
able identify tidal disruption events in the nuclei of galaxies, and trigger multiwavelength follow- 
up imaging and spectroscopy. 



Strubbe & Quataert (2009) [|48|] calculate models of tidal disruption flares at optical wavelengths 
in order to predict photometric and spectroscopic diagnostics. In these models, about half of the 
stellar debris accretes in a thick disk that radiates as a multicolor blackbody peaking around 0. 1 
keV. For ~ 1 month to ~ 1 year, the accretion rate remains high enough for the disk to shine 
at close to the Eddington luminosity. The other half of the stellar debris escapes from the BH, 
spreading into an angle that is larger for lower BH masses. This material is irradiated by the disk, 
and re-emits at longer wavelengths via recombination, bremsstrahlung, and lines. For BHs having 
Mbh ~ 1O^M0, the optical emission is dominated by the Rayleigh- Jeans tail of the disk's emission, 
and is very blue. For Mbh ~ 10^ — 10^ Mq, reprocessing by the unbound material can boost and 
redden the optical emission significantly. Super-Eddington rates of infall may also drive a wind 
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that could boost and redden the optical emission even more. 

In Figure 1, we use these emission models to predict detection rates with upcoming surveys. We 
assume a rate of 10^^ yr^^ per BH based on ROSAT detections [37] and conservative theoretical 
estimates, and assume it is independent of BH mass for simplicity. We use the BH mass function of 
Tundo et al. (2007) 1491], and assume it remains at lO^^Mpc"'^ down to 10^ Mq. Our calculations in- 
clude cosmological effects. We plot rates in the g and i band for the Pan-STARRS 1 Medium Deep 
Survey (PS 1 MDS) and 3n Survey, and for LSST. Detection rates peak at Mbh ~ few x IO'^Mq, but 
include a significant number of intermediate mass black holes (IMBHs) as well, particularly at i- 
band. With these assumptions, the MDS will detect 0.4 yr^^, Svr will detect 13 yr^^, and LSST will 
detect 130 yr^ This is in agreement with the predictions for Pan-STARRS and LSST in Gezari 
et al. (2008) [42] when we use the black-hole-mass dependent TDE rate used in their calculation 
BOll . Dust obscuration of the central regions of some galaxies might hide an appreciable fraction 
of events. Therefore synoptic searches in the infrared might be a useful complement to optical and 
X-ray studies. SASIR is a pre-phase A concept for a wide-field synoptic infrared survey with sci- 
ence operations expected at the end of the decade. Typical visits should yield depths of few /xJy in 
simultaneous YJHK imaging over 450 square degrees per night, with several thousands of square 
degrees that could be imaged on timescales relevant for TDE discovery. 

While surveys like Pan-STARRS 1 will confirm the feasibility of using optical surveys alone to 
discover and identify TDEs, wider, deeper surveys like LSST will be the factory that will produce 
large enough samples of these events to study their ensemble properties, and map out their rates as 
a function of galaxy type and redshift. TDEs can even be a powerful tool for studying the structure 
of galactic nuclei. The rate of TDEs is very sensitive to a galaxy's nuclear stellar density profile 
[30, 50] and non-axisymmetric geometry [la], as well as the presence of a nuclear star cluster [jsH]. 
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Fig 1. Predicted number of tidal disruption flare 
detections per year as a fimction of central BH 
mass for the Pan-STARRS Medium Deep (MDS), 
all-sky survey (Sir), and LSST survey. Adapted 
from Strubbe and Quataert (2009). 
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3.1 The Next Decade 



A more sophisticated theoretical treatment of TDEs than the invocation of simple fallback is clearly 
warranted. Just what are the signatures at radio and X-ray wavebands is clearly not well understood 
ab initio — and indeed the constraints placed across the electromagnetic spectrum will directly feed 
into and inform models of accretion in MBHs. 

It is clear that spectroscopic and multi-wavelength follow-up of optical detections will be im- 
portant. Optical/UV spectroscopy and color evolution will help confirm the nature of an event, and 
constrain the BH mass and pericenter distance of the disrupted star. These quantities are further 
indicated by the duration and bolometric luminosity of the flare, whose spectrum should peak at 
ultraviolet to soft X-ray wavelengths. By analogy to X-ray binaries, tidal flares may become bright 
at radio wavelengths at late times when the fallback rate subsides below ~ lO^^AfEdd, perhaps 
driving a jet [52]. The characteristic radius of the emission derived from the spectral energy distri- 
bution (SED) of a flare may be used to probe the innermost stable circular orbit (risco), which is 
strongly dependent on BH spin. 




Fig 2. Illustration of the observability of TDEs ( near 
peak brightness) to high redshift using JWST and TMT. 
The strong negative correction (K k, —7.5 logroll + z)) 
largely counteracts the di dimming and dust obscura- 
tion. 



Figure 2 illustrates how measurements of the broadband SED and detailed light curves of events 
could followed in the optical/NIR with instruments such as JWST and TMT, which will be sensitive 
to the Ray leigh- Jeans tail of the emission out to high redshifts. Indirect spectroscopic signatures of 
a TDE in the form of a "light echo" of fading high-ionization emission lines from the surrounding 
interstellar medium could be used as a diagnostic tool to study the energetics of the flare [47], 
while NIR imaging could be used to probe dusty gas in the cores of non-active galaxies that has 
been heated to high temperatures. 

The proposed EXIST mission plans to conduct a very sensitive all-sky hard X-ray time domain 
survey with a wide-field High Energy Telescope (HET) imaging the full sky every 3h in the 5 — 600 
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keV band but with TDE detections likely in the ~ 5 — 20 keV band. EXIST would have the unique 
ability to follow up all events with an onboard optical-Infrared Telescope (IRT) and Soft X-ray 
Imager (SXI). The discovery rate by HET of TDEs in the hard X-rays depends on the fraction 
of photons that are emitted at energies > 5 keV. However, evidence for a hardening with time of 



the spectra of two ROSAT TDE sources (RXJ1242.6 and NGC 5905) to T - 2.5 [j53|,l54l] points 
favorably towards the detection capabilities of EXIST iQ. Given that 30% of the ROSAT and 
XMM-Newton TDE candidates had F < 3 during the peak of the flare, and assuming that this is 
representative of the whole population, the BH mass function and tidal disruption rates described 
earlier would imply ~ 100 detections per year by the HET. 

A prompt (10 sec) flash of hard X-ray emission has been proposed to occur as a result of the 



tidal compression of the star at the time of disruption [15611 . Detection of this X-ray flash by the 
EXIST HET would mark the time of disruption to high accuracy, and the subsequent emission 
from the fallback of the debris could be monitored with IRT and SXI. Such precise measurements 
of the time of disruption and the shape of the power-law decay would place strict limits on the mass 
of the central BH and the mass of the star disrupted 

The proposed Wide-Field X-ray Telescope (WFXT) will survey the sky in the soft X-ray band 
(0.4 — 6 keV). While not explicitly a time domain survey, the mission plan includes a deep survey 
over 100 deg^ that will require multiple 1.5 ks visits to accumulate their deep exposure goal of 395 
ks. This would enable the discovery of tens of TDEs in an energy band that is expected to be close 
to the peak of their SEDs. 

We expect that by 2020, the systemic study of TDEs will be a major industry. A large 
census of events will allow us to map the masses and spins of SMBHs in galaxies as a function 
of redshift and galaxy type, as well as use the measured rates to probe the stellar structure of the 
host galaxy nuclei. This will allow the construction of an A^bh — cr* relation entirely independent 
from methods currently available, allowing for unique constraints on the central questions outlined 
concerning BH and galaxy growth. Although the detection sensitivity of the optical synoptic sur- 
veys shown in Figure 1 are "tuned" to BH masses of a limited mass range, (1 — 3) x IO^Mq, they 
will be a sensitive measure for deviations from the locally established scaling relations between the 
mass of the BH and their host galaxy properties. The low-mass tail of the detection sensitivity of 



LSST may even reveal the presence of IMBHs in dwarf galaxies and globular clusters pVQ. There 
is also the exciting potential to use the electromagnetic detection of a TDE to localize the burst 
of gravitational waves detectable by LISA when a star with a compact core approaches the tidal 
disruption radius or during the partial disruption of a white dwarf fl34l 56 . [58l1 . and to detect and 
localize recoiling SMBHs by TDEs from the cluster of stars which remains bound to them n5% . 
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